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I discuss recent developments in finite temperature lattice QCD, including the calcula- 
tion of the transition temperature, equation of state, color screening and meson spectral 
functions. 

1. Introduction 

In recent years considerable progress has been made in understanding QCD at finite 
temperature using lattice simulations. Calculations of bulk thermodynamic quantities 
have been done in the realistic case of 2+1 flavors with sufficiently small quark masses 
using different versions of improved staggered actions [ [TJ El El IH E] • This reduces or can 
even eliminate the uncertainty of extrapolations in the quark mass in different thermody- 
namics observables. 

One of the most prominent features of the deconfined phase of QCD is color screening. 
In non-Abelian gauge theories it is a very complicated phenomenon. Recently some 
progress in understanding color screening has been made by calculating the free energy 
of static charges on the lattice [El El E]- 

For many years lattice QCD did not provide any information on the properties of real 
time excitations at finite temperature. The situation has changed since few years and 
first calculations of meson spectral functions at finite temperature appeared [ El HH El 
El El El- In this contribution I am going to review some of these recent developments. 



2. Transition temperature and the equation of state 

One of the most interesting question for the lattice is the question about the nature 
of the finite temperature transition and the value of the temperature T c where it takes 
place. For very heavy quarks we have a 1st order deconfining transition. In the case 
of QCD with three degenerate flavors of quarks we expect a 1st order chiral transition 
for sufficiently small quark masses. In other cases there is no true phase transition but 
just a rapid crossover. Lattice simulations of 3 flavor QCD with improved staggered 
quarks (p4) using N t = 4 lattices indicate that the transition is first order only for very 
small quark masses, corresponding to pseudo-scalar meson masses of about 60 MeV [ 
IT7j . A recent study of the transition using effective models of QCD resulted in a similar 
estimate for the boundary in the quark mass plane, where the transition is 1st order [ 
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HH| . This makes it unlikely that for the interesting case of one heavier strange quark and 
two light u, d quarks corresponding to 140 MeV pion the transition is 1st order. However, 
calculations with unimproved staggered quarks suggest that the transition is 1st order for 
pseudo-scalar meson mass of about 300 MeV Thus the effect of the improvement 

is significant and we may expect that the improvement of flavor symmetry, which is 
broken in the staggered formulation, is very important. But even when using improved 
staggered fermions it is necessary to do the calculations at several lattice spacings in 
order to establish the continuum limit. Recently, extensive calculations have been done 
to clarify the nature of the transition in the 2+1 flavor QCD for physical quark masses 
using N t = 4, 6, 8 and 10 lattices. These calculations were done using the so-called stout 
improved staggered fermion formulations which is even superior to other more commonly 
used improved staggered actions (p4, asqtad) in terms of improvement of flavor symmetry. 
The result of this study was that the transition is not a true phase transition but only a 
rapid crossover [120]. 

Even-though there is no true phase transition in QCD thermodynamic observables 
change rapidly in a small temperature interval and the value of the transition temperature 
plays an important role. The flavor and quark mass dependence of many thermodynamic 
quantities is largely determined by the flavor and quark mass dependence of T c . For 
example, the pressure normalized by its ideal gas value for pure gauge theory, 2 flavor, 
2+1 flavor and 3 flavor QCD shows almost universal behavior as function of T/T c [I16j. 

The current status of lattice calculations of the transition temperature is summarized 
in Fig. la. Early calculations with pA action have been done for 2- and 3- flavor QCD 
for pion masses m v > 400 MeV [ [H] . Calculations with Wilson fermions have large 
uncertainties from chiral extrapolation. This is because the actual simulations are done 
in the region of large quark masses, corresponding to pion masses m, > 700MeV [ I2*2ll23j . 
There is also a calculation of T c using standard staggered fermion for the physical pion 
mass [123]. 

To better understand the systematics of the determination of T c in Fig. lb I show the 
transition temperature in units of the Sommer scale r as the function of the pion mass 
based on calculations on lattices with temporal extent N t = 4 and N t = 6 with p4 action 
[EJ. The pion masses in these calculations are smaller than 500 MeV with the smallest 
pion mass being around the physical value. Note that the value of T c calculated at two 
different lattice spacings are clearly different. The thin errorbars in Fig. lb represent the 
error in the determination of the lattice spacing a, i.e. the error in r /a. There is also an 
error in the determination of the gauge coupling constant (3 C = Q/g 2 . The combined error 
is shown in Fig. lb as a thick errorbar. For Nt — 4 calculations the error is dominated by 
the error in lattice spacing, while for N t = 6 it is dominated by the error in (3 C . With the 
data on r T c a chiral and continuum extrapolation has been attempted using the most 
simple Ansatz rQT c (m n , N t ) = r^T^^f 1 + A(r m n ) d + B/N 2 . From this extrapolation on 
gets the continuum value T c r = 0.457(7) [+8] [—2] for the physical pion mass m^r^ = 0.321 
[EJ. The central value was obtained using d = 1.08 expected from 0(4) scaling. To test 
the sensitivity to the chiral extrapolations d = 2 and 1 have also been used. The resulting 
uncertainty is shown as second and third error in square brackets. Using the best know 
value of r = 0.469(7) we obtain T c = 192(7) (4)MeV which is higher than the most of 
the previous values. It is also significantly higher than the chemical freezout temperature 
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Figure 1. Lattice data on the transition temperature from Refs. [ |2H d 13 13 1221 E3 and 
the chemical freezout temperature at RHIC [ 121)] (left). The transition temperature in 
units of the r T c from Ref. as function of the pion mass. 
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Figure 2. The energy density calculated with p4 and asqtad actions (left) and the pressure 
calculated with stout action. Calculation of the pressure with standard staggered action 
is also shown. 



at RHIC [|2H] which is shown in Fig. la as the vertical band. 

The equation of state has been calculated using both improved action and standard 
staggered fermion actions for different quark masses. The first calculation of the equation 
of state with improved action (p4) was done for Nt = 4 and at quite heavy quark masses, 
corresponding to pion mass of about 700 MeV at the transition temperature [I2E1- The 
quark mass was kept fixed in units of the temperature and not in physical units [ 26 . 
These may have an effect on the high temperature limit of the pressure and energy density 
[!2~71. More recent calculations with asqtad and stout improved action were done at smaller 
quark masses, namely for pions of about 200 and 300 MeV for asqtad and for physical 
pion mass for stout action. The result of these calculations is shown in Fig. 2. In Fig. 2a 
I show the energy density as function of T/T c for p4 and asqtad action. Calculations done 
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at Nt = 4 and N t = 6 agree well with each other. One can also see that the energy density 
calculated with the two different action agrees well despite the fact that the quark masses 
are significantly different. Close to the transition we may argue that the mass dependence 
of the energy density is taken care of by the mass dependnce of T c as discussed above. 

The pressure calculated with the stout action is shown in Fig 2b. Unlike p4 and asqtad 
the stout action impoves the flavor symmetry breaking but not the breaking of rotational 
symmetry which is important in the high temperature limit. This is probably the reason 
why the N t — 4 and N t = 6 results differ significantly for temperatures T > 1.3T C . In fact 
the cut-off dependence of the pressure for stout action in the high temperature phase seems 
to be as large as the cut-off dependence in the case of standard (unimproved) staggered 
action [EH- 111 Ref. [Hj it has been argued that the cutoff dependence can be taken care 
of by normalizing the pressure by the ideal gass value calculated on the finite lattice p§g. 
Indeed, p/p l gg shows significantly reduced cut-off dependence. However, calculations in 
pure gauge theory show that normalizing by the lattice ideal gas value may overestimate 
the cut-off effects [:29j. Therefore future studies are needed to settle this issue. 

3. Static quark and temporal meson correlators 

To study screening of static charges in high temperature QCD it is customary to calcu- 
late the color singlet Fi, octet F 8 and averaged free energies of a static quark anti-quark 

pair [mm eh 

e -F l( r,T)+C = I( TrW /(f)iyt( )), (1) 
o 

e -F(r,T) + c = h Tr w(f)TiW\0)) = ^exp(-F 1 (r,T)) + jj exp(— F 8 (r, T)). (2) 

While the color averaged free energy is gauge invariant, the singlet and octet free energies 
need gauge fixing. Usually the Coulomb gauge is used [ El 13 El El E21 E3] ■ There are also 
gauge invariant definitions of the singlet free energy in terms of Wilson loop [El 03] and 
eigenfunctions of the covariant Laplacian [ESI- m the high temperature leading order 
perturbation theory all definition give the same result. Also in the zero temperature limit 
different definitions should agree [ES]- The normalization constant C should be fixed such 
that the free energy of static quark anti-quark pair coincides with the zero temperature 
potential [Ej. In Fig.3a I show the results of the color singlet free energy in Coulomb 
gauge in 3 flavor QCD [ 9J. As one can see at very short distances the free energy is 
temperature independent and agrees with the zero temperature potential. It approaches 
a constant, F^, at large separation for all temperatures. However, as temperature is 
increasing F^ decreases as well as the distance where the free energy effectively flattens 
off. This can be thought of as manifestation of color screening. Similar results have been 
obtained in pure gauge theory [ El E2] and two flavor QCD [ EHj with the obvious 
difference that in pure gauge theory is infinite, as the free energy rises linearly with 
separation below the deconfinement temperature [E]- One can also calculate the entropy 
and the internal energy of the two static quarks, = —dF^/dT, = F^ + TS^. In 
Fig. 3b I show the entropy of two separated static quarks as function of the temperature 
for 2 and 3 flavor QCD as well as pure gauge theory. There is a large peak in Soo at the 
transition temperature. At large temperature seems to approach a constant value. 
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Figure 3. The singlet free energy in 3 flavor QCD [05] at different temperaure (left) and 
the entropy of isolated quark and anti-quark in pure gauge theory, 3- and 2-flavor QCD [ 
millini (right). 

The internal energy is very close to zero at high temperatures. The decrease of at 
high temperature is due to the presence of the entropy contribution, —TSoo, which at 
asymptotically high temperatures is given by — ^Ta s {mr, /T) ~ g 3 T and thus dominates 
the free energy [ 36 . Here mo is the leading order Debye mass. The internal energy is 
zero at order 0(g 3 ) and receives contributions from higher order. The above relation can 
be used for a non-perturbative definition of the Debye mass. 

The above analysis, done for the singlet free energy, can also be done for the color 
averaged free energy in the same way. At very small distances we have Fi = — (4/3)a s /r 
and F 8 = +(l/6)a s /r and thus F = F\ + Tln9 The color averaged free energy 

approaches the same constant F^, thus this quantity is gauge independent. However, the 
temperature dependence of F is much stronger than of F\ and one should consider much 
smaller distances to normalize it properly. This makes the lattice study of this quantity 
more difficult [ |B] . This is the main reason why most of the recent lattice calculations 
concentrated on the singlet free energy. 

Based on the simple picture of color screening it has been argued by Matsui and Satz 
that quarkonia should disappear at temperatures slightly above the deconfinement tem- 
perature To verify this conjecture and gain information about the quarkonium 
properties above the deconfinement transition meson spectral functions have to be calcu- 
lated. On the lattice this can be done using the relation between Euclidean and real time 
meson correlators, which in terms of the spectral function can be written as : 

G(r, T) = r dua(u, T)K(r, u>, T), K(r, u>, T) = ^f )} . (3) 

jo smh(o;/ 21 ) 

Given the data on the Euclidean meson correlator G(r, T) the meson spectral function can 
be calculated using the Maximum Entropy Method [EH]- For charmonium this was done 
by using correlators calculated on isotropic lattices [El E| as well as anisotropic lattices 
[El El El- Different charmonium states are identified as peaks in the spectral functions. 
It was shown that the spectral functions for the ground state charmonia change little across 
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Figure 4. The charmonium spectral function in the pseudo-scalar (left) and scalar (right) 
channels. 



the deconfinement transition [ [TTJ [TJJ EH EES] • This suggests that IS" (J/V, T] c ) charmonia 
can survive till temperatures as high as 1.6T C . Spectral functions of excited IP (x c ) states 
show large changes already at 1.1T C [[HUE! EI] • This is shown in Fig. 4 where the spectral 
functions in the pseudo-scalar and scalar channels calculated on anisotropic lattice [ ITS] 
are shown. We see no change in the pseudo-scalar spectral function which corresponds 
to i] c state. On the other hand the scalar spectral function corresponding to \c state 
shows dramatic modifications. Note that the spectral functions in Fig. 4 are zero above 
some energy. This is due finite lattice spacing which limits the maximal energy u. This 
behavior of the spectral function can be seen already at the level of the correlation function 
which does not require invoking statistical tools like MEM [ E3 El ESj • The pseudo- 
scalar spectral function show almost no temperature dependence across the deconfinement 
transition temperature, while the scalar correlator shows large enhancement over the 
zero temperature correlator [ El EI] • These features of the correlators are not easy to 
accommodate in potential model with screening [EH]. One of the problems with the 
spectral functions obtained using MEM is that it shows peaks even at energies uo > 5 
GeV where the continuum should dominate. This could be an artifact of the analysis 
method or the finite lattice spacing. Also the peaks at low uj appear to be quite broad 
(see Fig. 4) and thus may contain contributions not only from the ground state but also 
higher excited states. Clearly the spectral functions calculated on lattice cannot provide 
yet quantitative information about the charmonia properties at finite temperatures. The 
calculations of the charmonium spectral functions discussed above were done in quenched 
approximation. It will be interesting to see how these conclusions are modified when 
the effect os dynamical light quarks is taken into account. Preliminary calculation of 
the charmonium spectral functions in 2 flavor QCD suggest that the effect of dynamical 
quarks does not change the above conculsions significantly [ I42j . 

There also some preliminary studies of bottomonium correlators and spectral functions 
[HOIEI]- An intriguing outcome of these studies is the enhancement the scalar correlators 
above deconfinement temperature, similar to one observed in charmonium case. The scalar 
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bottomonium correlator corresponds to Xb state which is similar in size to IS charmonia 
states and thus is expected to survive above the deconfinement temperature. 

Transport processes should appear as peaks in the spectral functions at very small en- 
ergies. The width of these peaks is related to transport coefficients. Euclidean correlators 
are largely insensitive to the detail of these peaks and to very good approximation are 
determined only by the area under the peaks [ 03J HI] • Therefore there is little hope of 
getting transport coefficients from MEM. Nonetheless the presence of a transport contri- 
bution in the vector correlator has been identified [ I45j . 

Meson correlators and spectral functions have been calculated also for light quarks [ 
FTTH H51 Ej. These studies, however, are less systematic than for the charmonium case. 
One intriguing outcome of these calculations is the suppression of the vector spectral 
function at smaller energies which would result in the suppression of the thermal dilepton 
and photon rates. Clearly much more work is needed to verify or rule out this possibility. 
The lattice data on the vector correlator on the other hand show only small (about 10% ) 
deviation from the case of freely propagating quark anti-quark pair and can provide strin- 
gent constraint [QUI on the vector spectral functions. It would be nice to confront these 
findings with recent resummed perturbative calculations of the vector spectral function [ 

HE]. 

At the end of this section let me finally mentions that the temporal correlators and the 
corresponding spectral functions are useful not only to study mesonic excitations in high 
temperature QCD but also any real time excitation, including quark and gluon quasi- 
particles. This, however, requires gauge fixing. Preliminary results of lattice calculations 
of quark and gluon spectral functions in Coulomb gauge have been reported in [ I49j . 
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